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We give an Arakelov-theoretic definition of a naive height on divisors of 
degree zero on a hyperelliptic curve over a number field, and show that 
this naive height has computably bounded difference from the Neron-Tate 
height of the corresponding point on the Jacobian, a key ingredient being 
a theorem of Faltings and Hriljac (Theorem [2|) comparing the Neron-Tate 
height on the Jacobian to the arithmetic intersection pairing on the curve. 
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To simplify the exposition, we restrict to the case where our ground field is 
Q. 

We then use this result to give a new algorithm to compute the finite 
set of points on a hyperelliptic Jacobian of Neron-Tate height less than a 
given bound. This is important for a number of applications, for example 
to the problem of saturation (see |Sik95] o r [Sto 02]), to the computation of 
integral points on hyperelliptic curves (see [BMS "''08 ). to the use of Manin's 
algorithm |Man71] . and for numerically testing cases of the Conjecture of 
Birch and Swinnerton-Dyer. 

This paper bears some resemblance to the final two chapters of the au- 
thor's PhD thesis [Hoi 12b] . The author would like to thank Samir Siksek 
for introducing him to the problem, and also Steffen Miiller and Ariyan Ja- 
vanpeykar for some extremely helpful discussions, as well as very thorough 
readings of a draft version. 



1. Previous explicit computational work on Neron-Tate 

HEIGHTS 

The first definition of the Neron-Tate height was given by Neron in 1965 
|Ner65j . The problems of computing the height and computing sets of points 
of bounded height have been studied since the work of Tate in the 1960s (un- 
published), who gave a different definition from Neron which is sometimes 
easier to work with. Using this definition, Tate (unpublished), Dem'janenko 
|Dem68j . Zimmer |Zim76| . Silverman [Sil90j and more recently Cremona, 
Prickett and Siksek |CPS06j and Uchida [Uch06] have given increasingly re- 
fined algorithms in the case of elliptic curves. Meanwhile, in the direction of 
increasing genus, Flynn and Smart |FS97| gave an algorithm for the above 
problems in genus 2 building on work of Flynn |Fly93| , which was later 
modified by Stoll ([Sto99j and |Sto02| ). Recently, Stoll has announced an 
extension to genus 3 |Stol2j. 

The technique used by all these authors was to work with a projective 
embedding of the Jacobian or a quotient (usually the Kummer variety), to- 
gether with equations for the duplication maps, and thereby obtain results 
on heights using Tate's telescoping trick. However, such projective embed- 
dings become extremely hard to compute as the genus grows - for example, 
the Kummer variety is in genus 1, is a quartic hypersurface in for 
genus 2 and in genus 3 is given by a system of one quadric and 34 quartics 
m P^ [MuelO] . As such, it appears that to extend to much higher genus 
using these techniques will be impractical. 

In (Holl2a] . the author used techniques from Arakelov theory to give an 
algorithm to compute the Neron-Tate height of a point on a hyperelliptic 
curve, and a similar (though different) algorithm for the same problem was 
given by Miiller in |Muell) . Both gave computation examples in much 
higher genera (9 and 10 repectively) than had been possible with previous 
techniques, demonstrating that Arakelov theory can be effectively applied 
to this situation. This paper is a natural continuation of these ideas. 
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2. Other algorithms for heights in arbitrary genus 

As well as the development of practical algorithms to compute with 
heights on Jacobians of curves of low genus as described above, there is 
a parallel story of attempting to give algorithms which are valid for all (or 
large classes of) curves, though they may not be practical to implement. As 
with the computational approach, one usually starts by fixing a projective 
embedding of the Jacobian by a multiple of the theta divisor, and defines the 
naive height of a point to be the projective height of the image of that point 
under the embedding. Combined with bounds on the difference between the 
naive and Neron-Tate heights, this could be used to give an algorithm to 
compute the set of points of Neron-Tate height up to a given bound, and 
this is the approach that has been adopted by most authors mentioned in 
the previous section. 

It appears that it would be possible to give algorithms for bounding the 
difference between the Neron-Tate and naive heights for curves of arbitrary 
genus using this setup, but to the author's knowledge such algorithms have 
never been written down. Various authors have worked in this direction, and 
we mention certain of their contributions before comparing these results to 
those obtained in this paper using Arakelov theory. We restrict our attention 
to curves of genus at least 3, since the lower genus cases are covered above. 

Firstly, the works of Mumford |Mum66j and Zarhin and Manin |ZM72j 
describe the structure of the equations for abelian varieties embedded in pro- 
jective space and the corresponding heights and height differences, respec- 
tively. Several papers have been written on giving algorithms to construct 
'Mumford style' projective embeddings of Jacobians for curves of arbitrary 
genus. Most of these are restricted to the hyperelliptic case, for example 
the paper |VW98j of van Wamelen, and the PhD thesis [Rei72j of Reid, 
who consider embeddings with respect to 4??. The paper |And02j of Ander- 
son considers non-hyperelliptic curves, but only considers (slightly modified) 
embeddings with respect to 4'!9. Van Wamelen explicitly presents an algo- 
rithm, but the works of Reid and Anderson appear amenable to creating an 
algorithm. 

The other component needed is a bound on the difference between the 
Neron-Tate height and the naive height arising from such an embedding. 
Such a bound is given by Propositon 9.3 (page 665) in the paper |DP02j 
of David and Philippon, but requires that the Jacobian be embedded using 
16^9. 

The first consequence of this is that we cannot apply the results given 
above constructing projective embeddings, and so (to the author's knowl- 
edge) it is true to say that no algorithm for computing the set of rational 
points up to given Neron-Tate height has previously been given, even in the 
hyperelliptic case. On the other hand, it seems very likely that it would be 
possible with sufficient work to extend for example the work of Anderson to 
embeddings with respect to 16-d, though the result would be unlikely to be 
straightforward. 

The second consequence of the use of an embedding by 16?? in |DP02] is 
that the resulting algorithm is likely to be difficult to use in practise. Recall 
that the final step in finding the set of points of bounded Neron-Tate height 
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is always to search some region of projective space for rational points of 
naiVe height less than some other (related) bound. Using an embedding with 
respect to IGi? requires searching for points in projective space of dimension 
16^ — 1 where g is the genus of the curve, and so to search for points of naive 
height up to B will require time B^^^~^. 

2.1. Comparison of the algorithms. As discussed above, whilst (to the 
author's knowledge) no algorithm has previously been given to compute the 
number of points of Neron-Tate height up to a given bound, it is possible 
to make some remarks comparing the efficiency of a hypothetical algorithm 
based on projective embeddings with the algorithm given in this paper. We 
will thus assume for now that such a hypothetical 'projective embedding - 
based' algorithm is given, extending |And02] to give an embedding by 16??. 
Since the algorithm in this paper is only for hyperelliptic curves, we will 
necessarily restrict also to that case. 

2.1.1. Time taken to obtain a bound. It should then be possible to bound 
the time required to compute a bound on the difference between the nai've 
and Neron-Tate heights using a projective embedding, though the time re- 
quired would be likely to be large (consider for example that in practise it 
has not yet been possible to determine the equations for a projective embed- 
ding of the Jacobian of any curve of genus 3, or Kummer variety in genus 
4). 

In contrast, the algorithm presented in this paper uses at several points 
algorithms for which it is not possible at present to give bounds on the 
required run-time. In most cases this should not be too hard to remedy, 
but the step in Section [TOl which finds lower bounds on the values of theta 
functions on certain compact sets looks extremely hard to control. 

2.1.2. The size of the bound obtained. For similar reasons to those given 
above, it seems very hard to predict the size of the bound given by the 
algorithm in this paper. For a projective embedding, the bound is given in 
an explicit way in terms of the height of the image of the zero of the abelian 
variety under the projective embedding. It would be extremely interesting 
to know if this could be bounded a-priori. 

2.1.3. The search region. Recalling that the final stage of either algorithm 
is to search for rational points up to a given height in a projective space, 
we observe that the 'projective embedding' approach using 16-!? requires a 
search in a projective space of dimension 16^ — 1, whereas the method in this 
paper requires a search in projective space of dimension g. The low-genus 
cases suggest that this search may be the bottleneck of the algorithm, in 
which case this improvement is likely to be substantial in practise. 

2.2. Extension to the non-hyperelliptic case. It is reasonable to ask 
whether the algorithm given in this paper could be extended to the case 
of non-hyperelliptic curves. Sections [3lfT0] go through with relatively little 
change, the main difference being that a Merkl atlas should be constructed 
by pulling back along a Belyi map, as in [Javl2j . Section 13 is then redun- 
dant. Sections [11] and [12] do not carry through easily to the non-hyperelliptic 
case. It seems reasonable that it would be possible to obtain similar results 
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to those in these sections using naive heights on the image of the curve un- 
der the canonical embedding, but no details have been worked out by the 
author. 

3. Notation and choice of a resolution 

Throughout this paper, unless otherwise specified, C will be a hyperellip- 
tic curve of genus g over a number field K living inside weighted projective 
space P(l, 1, 5 + 1) with coordinates x, s, y. In Section [TT] we will need to 
specialise to -fC = Q, but until then we can work in full generality. We write 
Ok for the ring of integers in K. We assume C is defined by = f{x,s) 
where / = X]^£]'^ has integral coefficients /j. We write X = x/s, 

Y = y/s9+^, S = s/x and Y' = y/x9+^. 

Definition 1. For a number field L, a 'proper set of absolute values for L ' 
is a non-empty set of non-trivial absolute values on L such that the product 
formula holds. Note that we need absolute values (we do not allow their 
squares), and so if L is not totally real then a proper 'set' of absolute values 
is in fact not a set but a multi-set; we will ignore this distinction. 

These conditions determine a unique proper set of absolute values for each 
number field L, and we will denote it Ml. This uniqueness implies that if 
F/L is a finite extension, | — |^ an absolute value on L and \ — \^^ ; • • • ; l~L„ 
the absolute values on F extending v, then for all x ^ F we have (writing 
Np/L /o'" the norm from F to L): 

n 
i=l 

and for all x ^ L that 

(2) \Np,L{x)l = \{x)\f-'^^. 

We fix once and for all the following notation: 

• Mk a proper set of absolute values of K; 

• M% the subset of non- Archimedean absolute values of K\ 

• the subset of Archimedean absolute values of K; 

• K{iy) the residue field at an absolute value u G M^^. 

3.1. A regular model. Let 'i^'i denote the Zariski closure of C in Fc)j^(l, 1, 5+ 
1). A result of Hironaka, contained in his appendix to [CG084j (pages 102 
and 105) gives us an algorithm to resolve the singularities of by a sequence 
of blowups at closed points and along smooth curves (the latter replacing 
normalisations); we observe that 'Wi may locally be embedded in F^^^, and 
so Hironaka's result can be applied. 

We fix once and for all a choice of resolution of using the algorithm of 
Hironaka - thus we fix both the model ^ and the sequence of blowups used to 
obtain it. We will not need to assume that 'if is the minimal desingularisation 
of We also choose a component of the special fibre over each non- 
Archimedean absolute value, in order to fix the normalisation of the function 
$ later on. This is of course trivial over absolute values over which the 
special fibre is integral (for example, good reduction), and in practice will 
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be supplied by the 'rational point at infinity' on the odd-degree curves which 
we will restrict to in the final section. 

We write ly the usual (rational-valued) intersection pairing between divi- 
sors over 1/ G M% computed on (see |Lan88[ IV,§1]). 

Given an effective prime divisor D on C, we extend it to a horizontal 
divisor on ^ by taking the flat closure. We extend this to arbitrary divisors 
on C by linearity, and write D for this extension. 

4. A THEOREM OF FALTINGS - HrILJAC 

We shall make use of the following result which can be found in Lang's 
book |Lan88[ IV, §2]. 

Theorem 2. (Faltings |Fal84j and i7nl7ac [Hri83j ) Let D he a degree zero 
divisor on C , and let E he any divisor linearly equivalent to D hut with 
disjoint support. Then the height with respect to the "d-divisor of the point 
on Jac(C) corresponding to D is given hy 

h^{[0{D)]) = - J2 log\K{u)\L,{D + '^{D),E)-^ 9dAE) 

where ^ and go,u defined as follows: 

- $ sends a divisor of degree on the curve C to an element of the group 
of fihral Q-divisors on 'rf with order zero along the irreducihle component 
chosen ahove, such that for any degree-zero divisor D on C (with Zariski 
closure D) and fihral divisor Y on , we have ly [D + ^{D), = 0. 

- gD,u denotes a Green function with respect to the canonical (Arakelov) 
(1, l)-form for the divisor D, when C is viewed as a complex manifold via the 
emhedding v (see Section [9A\ for our normalisation of a Green's function). 

5. Outline 

Our definition of the height is analogous to the definition of the height of 
an element x in a number field K as 

(3) h(x) = ^ log+ \x\~^ . 

ueMK 

For each absolute value z/ of our number field, we will construct a metric or 
pseudo-metric di, on divisors which measures how far apart they are in the 
z^-adic topology. We then define 

(4) ^m)= E ^ogdyiD,Dr' 

ueMK 

where D' is a specified divisor which is linearly equivalent to —D. Since our 
curve is compact and our metrics continuous, the function d^^D, D')^^ is 
bounded below uniformly in D, and so we may use log in place of log"^. 

After setting up these metrics at non- Archimedean absolute values, we will 
first show how to bound the difference between the distance between two di- 
visors and their local arithmetic intersection pairing at a non- Archimedean 
absolute value. This will involve bounding the correction term <I> (see |Lan881 
III§3]), and then comparing distances between divisors and lengths of mod- 
ules, in Section [8l The hardest aspect of this will be allowing for the fact that 
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the model of C obtained by taking the closure inside projective space over 
Ok is not in general a regular scheme, so we must compute precisely how 
the process of resolving its singularities will affect the intersection pairing. 

We then use a theorem of Merkl to construct a weak-pseudo-metric on C 
at each Archimedean absolute value, and to bound the difference between 
this fTinction and the local Ncron pairing. Finally, wc write down two more 
nai'vc heights, with successively simpler definitions, each time bounding the 
difference from the Neron-Tate height. We give an algorithm to compute the 
number of points of bounded height for the simplest of these naive heights, 
completing the algorithm. 



6. Metrics on C 
We begin by setting up a collection of metrics. 



Definition 3. For each absolute value u G Mk, we define [K^^ ,\ — \) to 
he an algebraic closure of the completion Ky together with the norm which 
restricts to v on K C Ki^^. For non- Archimedean absolute values v we 
define 

(5) d, : C{K^^3) X C{Ki^9) ^ 

by 

d;/((a;p : Sp '. yp), {Xq '. Sq '. yq)) 



= < 



max ( \Xp/ Sp Xql Sq\ , 
max ( \Sp/Xp — Sq/Xq \ , 



yp/sp'^ - yqls%^^ \ if \xp\ < \sp\ and \xq\ < 

Vp/xp^^ - yq/xp^ ^ if \xp\ > \sp\ and \xq\ > 

otherwise. 



Proposition 4. For each v G M^^, d = is a metric on C{Kf'^). More- 
over, for each such v, we have di^{p, q) < 1 for all p and q. 

Proof. Only the triangle inequality is non-obvious. We begin by observing 
that ii{x:s:y)e C{K^^^) then 



(6) |x| < |s| 
and 

(7) \x\ > \s\ 
The first implication holds as 



\y\ < \s\'^' 

\y\ < 



(8) 



\y\ 



^ fiX'S 



2g+2-i 



< ( max|/i 



and the other case is similar. Combining this with the fact that | — | is p-adic, 
we see for all p, q E C{K^ ^) that d{p, q) < 1. 

From now on we proceed case-by-case. Let p = {xp, Sp, yp), q = {xq, Sq, yq) 
and r = {xr,Sr,yr)- Note that the change of coordinates x i-)- s, s i-)- 
X, y/s^^^ I— )• y/x^~^^ preserves the metric (replacing / by its reciprocal 
polynomial). As such, we may assume without loss of generality that |xp| < 
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Case 1: \xq\ < \sq\ and \xr\ < \sr\- Then 



(9) 

d{p,q) +d{q,r) 



max 



> max 



> d{p, r) 



Xp 


Xq 






Sp 


Sq 




P 


Xp 
Sp 


Xq 
Sq 


+ 


Xq 

Sq 



^ + max ^ 


X q Xf 










Sq S^ 





Vq Vr 



S^- 



Vp Vq 



,9+1 ,9+1 



+ 



Vq Vr 



,9+1 ,9+1 



Case 2: \xq\ < \sq\ and \xr\ > \sr\- If \xq\ < \sq\ then 

(10) d(p, q)+d{q,r) = d{p,q) + l>l> d(p, r) , 

so we may assume \xq\ = \sq\. Now if = \sp\ then we are back to Case 1, 
so assume \xp\ < \sp\. Then since | — | is p-adic, we have \xp/sp — Xq/sq\ = 
max(|a;p/sp| , \xq/sq\) = 1, so 

(11) 

d{p,q) +d{q,r) 



> max ^ 


Xp 


Xq 


+ 


X q X y 






Sp 


Sq 




Sq Sr 





Vp Vq 



,9+1 ,9+1 



+ 



Vq Vr 



,9+1 ,9+1 



> 1 > d{p, r). 



Case 3: \xq\ > \sq\ and \xr\ < \sr\- If \xp\ < \sp\ then 

(12) d{p, q) + d{q, r) = 1 + d{q, r) > 1 > d(p, r), 
and if \xr\ < \sr\ then 

(13) d(p, q) + d{q, r) = d{p, q) + 1 > 1 > d{p, r). 

Otherwise we are back to Case 1. 

Case 4: \xq\ > \sq\ and \xr\ > \sr\- Interchanging p and r reduces us to 
the second case. □ 



7. NoN- Archimedean I: the $ term 

In this section we work locally over a non-Archimedean absolute value, 
so for the remainder of this section let denote a regular model of the 
curve C over a discrete valuation ring R finite over Zp. We replace R by 
an unramified extension such that all irreducible components of the special 
fibre of ^ over R are geometrically irreducible. If is smooth over R, then 
the constructions in this section are trivial. We will show how to bound the 
values of the function as defined in |Lan881 III§3]. 

Let F denote the free abelian group generated by prime divisors supported 
on the special fibre, and let V denote the finite-dimensional Q-vector space 
obtained by tensoring F over Z with Q. Let M : y x y — )■ Q be the map 
induced by tensoring the restriction of the intersection pairing on '^^ to its 
special fibre with Q. Then V has a canonical basis of fibral prime divisors, 
so we may confuse M with its matrix in this basis. Call the basis vectors 
Yi . . . Yn, we use the same labels for the corresponding fibral prime divisors. 
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We will make use of the Moore-Penrose pseudo-inverse (first defined in 
[Pen55j ) of the matrix M, which we denote . We do not need the full 
definition, only existence and the fact that if for any matrix A the linear 
system Ax = b has any solutions, then a solution is given by x = A^b. 

We will also need the definition of a reduced divisor on a hyperelliptic 
curve: 

Definition 5. We say that a divisor D on C is semi-reduced if it is effec- 
tive, any Weierstrass point appearing in its support has coefficient 1, and 
if its support does not contain any pairs of points p, ( where p~ denotes 
the image of p under the hyperelliptic involution). If in addition we have 
deg(-D) < g, then we say D is reduced. 

Proposition 6. Let Af + denote the Moore-Penrose pseudo-inverse of M , 
let ru- denote the infimum of the entries of and m+ their supremum. 
Let D = Di — D2 and E = Ei — E2 be differences of reduced divisors on 
C with no common points in their supports, and assume that D and E both 
have degree zero. Then 

(14) \L,{^{D),E)\<g\m+-m^). 

Proof. Let d denote the vector X^ILi '•i^ {D, 5^,) Yi, and similarly set e to equal 
Y17=i (-^'^j) ^) ^ P^i^ of vectors in V. Now by definition of <I> we have 
that for all vectors v G V: 

(15) V d'^ + V M ■ <S>{Df = 0, 
and hence that 

(16) d'^ = -M ■ <^{Df. 

Using the property of the Moore-Penrose pseudo-inverse given above, we 
can take ^{D) to be —d ■ (M"*")"^, and so we find 

(17) Lu{'^{D),E) = -d- {M+f -e^. 

Now since D and E are differences of reduced divisors, d and e are vectors 
each formed by assigning g copies of '+1' and g copies of ' — 1' to the basis 
elements Yi,. . . ,Yn (allowing multiple ±ls to be assigned to a single basis 
vector), and so the result easily follows. □ 

Definition 7. Using the above proposition, we can define a computable con- 
stant depending only on ^ such that for all differences of reduced divisors 
D and E with no common points in their supports, we have 



(18) 



^ i4^,{D),E) 



< 



»1. 



8. Non-Archimedean II: local comparison of metrics and 

INTERSECTION PAIRINGS 

We compare the metrics and intersection pairing for points on C. The 
main difficulty is that of working with regular models for C; the naive pro- 
jective closure of the generic fibre is not in general regular, but rather must 
be modified by a sequence of blowups at smooth centres to obtain a regular 
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model. We must determine how these modifications wih affect the inter- 
section numbers, and also keep careful track of the base field since regular 
models are not in general stable under ramified base change. 

Let b denote the longest length of a chain of blowups involved in obtaining 
'if from (one blowup is considered to follow another if the centre of one 
blowup is contained in the exceptional locus of the previous one). This is 
finite since Hironaka's algorithm terminates in a finite number of steps. Note 
that we fixed a choice of the resolution and the algorithm, so that and 
b are well-defined. 

Recalling the definition of the metric d(— , — ) from Section [6l the aim of 
this section is the following global result (which will follow easily from the 
subsequent local results): 

Proposition 8. There exists a computable constant SS2 such that for all 
reduced divisors D and E (with no common points in their supports) on C 
with closures Si and S on , we have 



where L/K is any finite extension over which both D and E have pointwise 
rational support. 

Proof. This is a global version of Proposition [21 whose validity we shall 
assume for now. We begin by noting that there are only finitely many 
absolute values of bad reduction for K, and that at any absolute value of 
K of good reduction, we have 



by applying Proposition [9] in the case bi, = 0. 

As such, it suffices to fix a absolute value u of K, and then bound the 
quantity 



(where the sums are over absolute values uj oi L dividing v) uniformly in D, 
E and L. We begin by noting that the expression does not depend on the 
choice of L. Now L ®k is a product of finite extensions of the completion 
Ky^ and moreover the degree and number of the extensions may be bounded 
in terms of the genus g of C . We may then apply Proposition [9] to each of 
these fields to obtain the result. □ 

The remainder of this section will be given over to stating and proving the 
local result Proposition [9j For this, we fix a non- Archimedean absolute value 
V of K. By base change, we have a chosen resolution '^k^ of the closure of 
C in weighted projective space over Ok^- Let by denote the longest length 
of a chain of blowups involved in obtaining this resolution (one blowup is 
considered to follow another if the centre of one blowup is contained in 



(19) 




(20) 




(21) 
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the exceptional locus of the previous one) This is finite since Hironaka's 
algorithm terminates in a finite number of steps. Note that b,^ = ii is 
smooth over Ok^- 

For the remainder of this section, let D and E be effective divisors on 
C with disjoint support, of degrees d and e respectively. Let L^/Ky be a 
finite extension of degree m with residue field I such that D and E are both 
pointwise rational over L^. Write D = Yli=iPi ^^'^ ^ — Yli=i 1i- Write to 
for the maximal ideal of Ol^- 

Proposition 9. Let S> and S denote the Zariski closures of D and E re- 
spectively on the regular model over Ok„ ■ Then 

(22) -b,de<i,{&,<^)- Vl ^°g^(rr^ vl^O- 

To avoid an excess of notation, we will from now on drop the subscript v 
from the fields and models we are considering - no confusion should result, 
since we will exclusively be working locally. 

To prove Proposition [9l we will need a sequence of lemmas: 

Lemma 10. Let p ^ q e C{L) = HomL(L, Cl). Write 

(23) length^,, (^^L^^^^) , 
f^l^ ^ + ^ 

where the sum is over closed points Vt of 'XOk lying over uo, and Lp 
and Lq are defining ideal sheaves for the closures in ^\ x Ol of the images 
of p and q in C L. Then 

(24) Lp^q = mlog 



d{p, q) 



Proof. Write p = {xp : Sp : Up), q = {xq : Sq : Vq) with Xp, Sp, Xq, Sq G Ol- 
If |xp| < \sp\ and \xq\ > \sq\ or vice versa, then p and q do not meet on the 
special fibre so ip,q) = 0, and by definition we see that dp{p,q) = 1. 

Otherwise, possibly after changing coordinates, we may assume that p 
and q are of the form (xp : 1 : ?/p) and (xg : 1 : ?/g) respectively, for Xp, i/p, 
Xg, yq G Ol- Writing F for the (integral) defining equation of C on the 
coordinate chart containing p and q, and taking to be the closed point 
where p and q meet, we have 

Ovixoj^OL,n ^ OL[x,y](^^^y) 



(25) 



so 



Ip~\~Iq {T, X Xp,y ypiX Xq^y yq) 



(Xp Xq,yp yq) 



(26) lengthc)^ ^_^o<o^Ol^^ ^ ^.^ (orda,(2;p - Xg), ord^(yp - yq)) . 

Now given a £ L, we find 

(27) log(#/)ord^(a) = -mlog|a| , 
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SO 

l^r^o-fh f ^'^iXoK(^L,n \ __ mm{-log\xp-Xq\,-log\yp-yq\) 
and hence 

(29) Ip,g = mmin(-log|xp - Xg\ , - log |yp - yq\) . 
However, 

(30) log(l/d(p,g)) = m.m{-log\xp - Xg\ , -log\yp - yg\) , 

so we are done. □ 

Lemma 11. Recalling that over L we can write D = Yli=iPi ^''^^ ^ ~ 
Yli=i 1i' define Uij to be the closed point of^i 'Xqk where pi meets 
Qj if such exists, and the unit ideal otherwise. Letting Id and Ie denote the 
ideal sheaves of the closures of D and E respectively on we have 

(3. E.-'V (^) -en.,., (^.i;^^^ 

The analogous statement on ^ also holds. 

Proof. We may decompose Id and Ie into iterated extensions of the sheaves 
/p. and Iq. , whereupon the result follows from additivity of lengths in exact 
sequences. □ 

Lemma 12. Let lo and Ie denote the ideal sheaves on corresponding 

to the closures of the divisors D and E respectively. 

(32) 

(^^) • = ^^"^^"^^ {{Io+lE)l2o,) ■ 

The analogous statement on ^ also holds. 

Proof. Let M be a finite length Ox-module. We show 

(33) lengtho^(M) • ram. deg(L/i^) = lengthy (S)Ok ^l)- 

Let M = Mq C Ml C • • • C = be a composition series for M, so 
each Mi/Mi^i is simple. Since Ok is local, we have by [MatSOl pl2] that 

(34) M,/Mi+i ^ Ok/xxik. 
By additivity of lengths, it suffices to show 

(35) lengthen ^i) = ram. deg(L/if), 

but this is clear since xuk ■ Ol = xn^^''^^^''^^\ □ 

Lemma 13. Let </> : — )• '^2 be one of the blowups involved in obtaining 
from Let p Cl{L). Then 

(36) < lengthc,^ ( V^^) '^''''^^^^^ i ^If + h ) ~ ^^^^S^^/^)- 
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Proof. In this proof, we will omit the subscripts 'Ol' from the lengths, since 
all lengths will be taken as O^-modules. If p does not meet g on ^2 x Ol 
then both the lengths are zero, so we are done. Otherwise, let be the 
closed point on ^2 x Ol where p meets q, and let a be the closed point of 
"^2 such that 0, lies over a. 

Let R denote the local ring of the (three-dimensional) ambient space to 
^2 at ct, and similarly let A be the local ring of ^2 at ct. Let B G Rhe the 
centre of the localisation of (p at a. After etale base-change, we may assume 
that we have 

(37) R = OK[[x,y]](x^y^a) 

where Ok is finite etale over Ok and a is a unif or miser in Ok, and that 

(38) B = {x,y,a) or B = {x,a), 

depending on whether we are blowing up a point or a smooth fibral curve. 

Blowups commute with flat base change, and the strict transform of a 
closed subscheme under a blowup is the corresponding blowup of that closed 
subscheme (see |Liu021 Corollary 8.1.17]), so we can be relaxed with our 
notation. We may write 

(39) p={x- axp, y - ay^) q = {x- axg, y - ayq) 

where Xp, yp, Xq and yq are in Ol ■ Ok- Setting uj' to be a uniformiser in 
the maximal ideal of Ok ■ Ol, we have 

(40) length ( ^^^7^ ) ^ "^^'^ (ord^/ (axp - axg), ord^/(ayp - ayq)) . 

In the case B = {x,y,a) we look at the affine patch of the blowup given by 
setting a 7^ 0; the equations for p and q transform into 

(41) p' = (x- Xp, y - yp) and q = {x - Xq, y - yq), 
so 

length ( ^^X7^ ) ^ min(orda;'(xp - Xq),OT:d^>{yp - yq)) 

(42) V + ^9 / 

= length -ord^,(a). 

In the case B = {x,a) we look again at the affine patch of the blowup given 
by setting a 7^ 0; the equations for p and q transform into 

(43) p' = (x- Xp, y - ayp) and q' = {x - Xq, y - ayq), 
so 

length ( ^^xf^ ) = min (ordo;' (xp - Xg), ord^^/ (ayp - ayq)) 

(44) \h + ^gJ 

= length (j^fj^) - (0 or 1) ord,,(a), 

so the result follows from the fact that, since Ok is unramified over Ok, we 
have 

(45) ordi^'(a) = ram. deg(L • K/K) = ram. deg(L/i^). 



14 



DAVID HOLMES 



□ 

Proof of Proposition To prove Proposition we apply Lemmata [TDl 
[IT] and [T2] in that order to find that there exists < /3 < bde such that 

(46) 



Ip + Ig 



= -log(#KM) lengthy, +/3 

1 f \ 

= — log(#K(w)) lengthy - — -— • ram. deg(L/i<") + /3 
m \Id + IeJ 

□ 

Remark 14. A'^oie i/iai /3 = for all hut finitely many absolute values v. 

9. Merkl's Theorem 

There are two key tools which we will need to handle local heights at 
Archimedean absolute values. The first is a theorem of Merkl, and the 
second is a formula of de Jong taken from his PhD thesis. In this section we 
will state Merkl's theorem |CE^11 , along with improvements due to Peter 
Bruin [Bru06] . 

9.1. Statement of the Theorem. Let X be a compact connected Rie- 
mann surface of positive genus and let ^ denote the canonical (Arakelov) 
(l,l)-form on X. We write g{—, —) for the Green's function on X x X as- 
sociated to and G : X X X — 7- M>o for its exponential, which we normalise 
to satisfy the following three properties: 

1) G{p, q) is a smooth function on XxX and vanishes only at the diagonal. 
For a fixed p G X, an open neighbourhood U oi p and a local coordinate z 
on U we can write 

(47) logG{p,q) =log|z(g)| +a{q) 

for p ^ q £ U , with a a smooth function 

2) for all p £ X we have dqdg log G{p, qY = 2Trifi{q) for q p. 

3) for all p G X, we have 

(48) / logG{p,q)fi{q)=0. 

Jx 

We give a trivially modified version of the definition of a Merkl atlas which 
will be more convenient for our purposes - the usual version has ri = r2, 
but the additional data we associate will be useful later on: 

Definition 15 (Merkl atlas). A Merkl atlas for X is a quadruple 

(49) ({C/„z,-)}^=i,ri,r2,M,ci), 
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where {{Uj,Zj)}^^^ is a finite atlas for X, and 1/2 < r2 < ri < 1, Af > 1 
and ci > are real numbers such that the following properties hold: 
1) Each Zj{Uj) is the open unit disc. 



2) The open sets C/J^ {x G Uj : \zj{x)\ < with 1 < j < n cover X. 

3) For all 1 < < n, the function \d Zj / d Zj'\ on Uj Pi Uj' is bounded 
from above by M . 

4 ) For I < j < n, write fj. = iFj d Zj Ad Zj on Uj . Then < Fj (z) < ci 
for all z £ Uj. 

Theorem 16. Let M = ({^7,, n, r2, M, ci) be a Merkl atlas for X. 

Set 

330n 1 

(50) =^(M) = — log + 13.2nci + (n - 1) log M. 

(1 — n)'^/"' 1 — ri 

Then 

(51) sup g{p,q) < jSg{m). 

{p,q)<=XxX\A 

Moreover, for every 1 < j < n and p ^ q G U^ , we have 

(52) \g{p,q) - \og\zj{p) - Zj{q)\\ < ^(M). 

Proof. The usual version of the Theorem has ri = r2, but the proof of this 
version follows trivially from the usual one. Merkl proved a less explicit 
version of this theorem in CE"'"ll[ Chapter 10]. A version similar to this 



was given by Peter Bruin in his Masters Thesis jBru061 Theorem 3.1], and 
the result in precisely this form is proven in Bruin's appendix to |Javl2] . □ 

9.2. Constructing a Merkl atlas. For the remainder of the section, we 
fix a complex absolute value of our ground field K, so that we may by slight 
abuse of notation view a smooth variety over K as a complex manifold. 

We will need to explicitly compute the constant =^(M) for our hyperellip- 
tic curve C. In order to do this, we must first fix a construction of a Merkl 
atlas for C, which we do using ideas from |Javl2) . Recall that we write 
TT : C — 7- for the projction map. We will construct an atlas on P^, then 
lift it to an atlas for C. 

9.2.1. An atlas for P^. We first construct a finite cover of the Riemann 
sphere as a metric space (with the chordal metric) by open discs with the 
property that for any two such disks Di and D2, if Di contains the centre 
of D2, then Di = D2. An algorithm to achieve this is given in Appendixll3l 
Next we need to relate this topological cover to a holomorphic cover. We 
do this one disk at a time; 



x,s 



Definition 17. given a disk D with centre p, translate p to {0 : 1) £ P], 
using a unitary transformation. Then set Z£){q) = x{q)/\x{qQ)\ for any qo 
lying on the boundary of D. We set R = \x{qQ)\ for any qo on the boundary 
of D . Set r2 = \qi\ where qi is any point in D lying at distance r from p in 
the chordal metric, and pick any r2 < ri < 1. 

Because of the way we set up this cover of P^, it should be clear that 
it satisfies the criteria (1), (2) and (3) in the definition of a Merkl atlas. 
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However, unsurprisingly our real interest lies in constructing a Merkl atlas 
for C, which we will achieve by lifting this atlas for P^. 

9.2.2. Lifting to a Merkl atlas for C . There are two cases to consider in 
lifting disks D from U to C. The first and easiest is when D is not centred at 
(and therfore does not contain) a branch point. Then the preimage ■k~^{D) 
is a pair of disjoint open sets C/i, U2 in C (since vr is etale outside i?), and 
we take these open sets together with the maps zu = zd ° obtained by 
composing the projection to with the map zd '■ U ^ C 

The harder case arises when D contains (and so is centered at) a branch 
point b. Writing Bq{1) for the unit disk in C, we take the unique map z 
making the following diagram commute: 



D ^-Bo(l) 

ZD ' 

9.3. Making the constants in Merkl's theorem explicit. Now that 
we have shown how to construct a Merkl atlas, we need to show how to 
explicitly compute the real numbers M and ci appearing in the statement 
of Merkl's theorem. 

9.3.1. The computation of M. It suffices to give an algorithm to find an 
upper bound on jdzj/dzjl on Di n Dj for fixed charts {Di,Zi) and {Dj,Zj) 
in the Merkl atlas. If neither Di nor Dj contains a Weierstrass point, write 
Pi and pj for their centres. Then 

(53) z,{X)-^~^^^ 



and so differentiating we find 



(54) ^ = ^^4^ 

so to find a bound on the absolute value of the derivative d Zj/ d Zj on DiCiDj, 
it suffices to bound this rational function on the domain 

(55) {X £C: \MiX\ < R and \MjX\ < R}, 
where 



XXp^ + 1 



XXp, + 1 



(56) Mi 



Pi 



is the unitary matrix translating pi to X = 0, and R is as in Definition [T71 
By our construction, at most one of Di and Dj can contain a Weierstrass 
point if Di n Dj is non-empty. In this case, we assume without loss of 
generality that Di contains the Weierstrass point, and we use a unitary 
transformation to move its centre pi to X = 0. We proceed in a simiilar 
manner to before, using the above descriptions of the matrices Mi and Mj. 
Zi and Zj are the coordinates on the charts Di and Dj given in Section [9.2.21 
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(so Zj is in effect the same coordinate as before, but Zi is the square-root of 
the coordinate Zi used in the previous case). Differentiating, we obtain 

(57) dzi _ 1 \Xpf + l {XX^-^+lf 1 



dzj 2|x^j2^i(XXp, + l)3/4X-X,/ 
whose absolute value we bound on the same domain as before. 

9.3.2. The computation of c\. Since we may translate any Merkl chart to 
be centred at X = using a unitary matrix, we will treat just this case. 
We may thus write our chart U = {{X,Y) £ C : \X\ < R} and z = X/R 
OY z = \/X/R, depending on whether U contains a Weierstrass point. We 
know that a basis of differential forms on C is given by 

, , X"dX 

(58) — - — < n < 5 - 1. 

Petersson inner products may be computed to any required precision by nu- 
merical integration (though care must be taken to either bound the deriva- 
tives, use interval arithmetic or similar in order to yield a provably correct 
result), and applying Gram-Schmidt so we can express the Arakelov (1,1) 
form in terms of this basis. It thus suffices to bound the functions -Fm,n 
defined by 



(59) — - — = F^,„dzAdz. 

We consider first the case where U does not contain a Weierstrass point. 
Then 

(60) dz^dz = ^AX ^dX, 

R 

so 

v"-*-) i'm.,n — \^Y^ ' 

which is bounded on U. 

If U contains (and so is centred at) a Weierstrass point, we have z = 
\JXIR, so 

(62) dzAdz = \ dX AdX = -^dX AdX, 

so 

41X1 X™X" 



(63) = 



R \Yf ■ 

Now since = F{X) = XF^iX), this reduces to 

4 



(64) F, 



m.n 



R\Fo{X)y 
which is bounded on U. 
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10. Theta functions to compute Green's functions 

The formulae and bounds given for Green's functions in the previous 
section only apply to pairs of points on C which are sufficiently close together 
- more precisely, to pairs of points which can be contained in the same chart 
of our Merkl atlas (more precisely, in U'^'^ where C/ is a chart). To handle 
the remaining points, we will use a formula from the thesis of Robin de Jong 
( |dJ04] and |DJ05j ). which expresses Green's functions in terms of theta 
functions. In order to give uniform bounds, we will use this formula first 
to evaluate the Green's function at a finite set of pairs of points, and then 
to bound its derivative so as to obtain bounds for points nearby. C will be 
a hyperelliptic curve of genus g (though de Jong's results apply to general 
Riemann surfaces). 

10.1. Pseudo- metrics to approximate Green's functions. In order 
to effectively use de Jong's results on theta functions, and to combine them 
with the bounds given by Merkl's theorem, we will need to define a weak- 
pseudo- metric d on points of C such that log(d(pi,p2)) is 'not too far from 
a Green's function'. This function d is analogous to the metrics djy given for 
non- Archimedean absolute values v, and it will be combined with them in 
the definition of the naive height. 

We wish to be able to use Merkl's Theorem to show that our weak-pseudo- 
metric is close to a Green's function for points that are 'close together', and 
also to be able to apply results with theta functions to show that our weak- 
pseudo- metric is close to a Green's function for points that are 'far apart'. 
To enable this, we begin with a rather contrived definition of a weak-pseudo- 
metric. 

Definition 18. Let M 6e a Merkl atlas for a hyperelliptic curve C con- 
structed as above and with r2 < ri. Let < p < {ri — r2)/4, and let V' be a 
finite cover of by closed disks of radius p in the chordal metric, with the 
property that if any disk contains a branch point then it is centred at that 
branch point. Lift V' to a cover V of C just as we did to obtain the Merkl 
atlas. 

Given two charts V\, V2 in V with centres qi, q2 respectively, define 
d{Vi,V2) = G{qi,q2) (where G is the exponential Green's function) if ViCl 
V2 = 0, and d{Vi, V2) = 00 otherwise. 

Given two pairs {pj , C/J^ ) with Uj a chart from M and pj £ Uj^ , define 

l^iiPi) - ziip2)\ ifUi = U2 
00 otherwise 

Given two triples {pj,Uj,Vj) with pj £ [7j^, pj G Vj, Uj a chart from the 
Merkl atlas M and Vj a chart from the cover V, we define 

(66) d{{pi,UuVi),p2,U2,V2)) = mm{d{{pi,Ui),{p2,U2)),d{VuV2) 
Given two points pi, p2 £ C , we define 

(67) d{pi,p2) = min d((pi, C/i, Fi), (p2, ^^2, ^2)), 



(65) d((pi,C/i),(p2,^72)) 



where the minimum is taken over all charts Ui, U2, Vi, V2 such that pj G 



Uj^ n Vj as before. 



HEIGHTS ON HYPERELLIPTIC JACOBIANS 



19 



It is clear that the function d will depend on the choice of p and also on 
the choice of the cover V. We will sometimes write dp or 'ddpy when it is 
important to make this explicit. 

Lemma 19. For all < p < (ri — r2)/4, the function dp : C x C — )• 
M>o U {oo} takes finite values. 

Proof. It suffices to show that for every pi, p2, we have either 

1) there exists a chart U gM. such that pi € C/*"^ and p2 G U^^. 

2) There exist disjoint Vi, V2 GV such that pi G Vi and P2 ^ V2. 

Suppose (2) fails. Since the V cover C, there must exist a pair of over- 
lapping charts V, V £ V such that pi £ V and p2 € V. In particular, 
in the chordal metric on we find d(7r(pi), -7r(p2)) < ?'i — i"2. Since the 
U"^^ cover C, there exists a U £ U such that pi G U^^. Since the chordal 
distance d(7r(pi), 7r(p2)) < — '"2, we see n{p2) £ -7r{U"^^). Now if [/ is a 
chart containing a Weierstrass point it is clear that p2 £ U'^'^ , and if not 
then since the pj are contained in overlapping charts V, we see again that 
P2£U''\ □ 

The aim of the remainder of this section is to prove the following result: 

Lemma 20. For all pi, P2 £ C, and for all < p < (ri — r2)/4, we have 

(68) \logdp{pi,p2)-g{pi,P2)\ <max(^(M),c(p)), 
where the constant c{p) is as defined in Proposition {2^ 

Definition 21. To ease the notation, we will from now on fix once and for 
all a Merkl Atlas, a constant < p < {ri — r2)/4 and cover V. 

We can then define a constant ^3 =^ max{^(M.),c{p)), for use later on. 
The idea of the proof is as follows: 

firstly, if we are in situation (1) of the proof of Lemma [T9l then the bound 
follows immediately from Merkl's Theorem; 

secondly, if we are in situation (2), then we need to bound the difference 
between the value of the Green's function at the centres of the charts Vj with 
the value of the Green's function evaluated at any points in those charts. To 
do this, we bound the derivative of the Green's function via a relationship 
with theta functions. 

10.2. Some results of de Jong. We state the main definitions and results 
we will need from jd J04] . We restrict to the hyper elliptic case and genus 
5 > 2 as this simplifies the statements somewhat, though the results hold 
more generally. We take a Weierstrass point 00 of the curve as a base-point 
for the Abel-Jacobi map a to the analytic Jacobian. 

Definition 22 ( |dJ04l Proposition 1.4.4]). Let r £ %g lie in the Siegel upper 
half space, and z £ . We then define 

(69) p||(z;r) = (det(9(T)))i/^ exp (-^9(z)*(9(r))-i9(z)) \^{z;t)\. 

Definition 23 ( |dJ041 Definition 2.1.1]). We define the invariant S{C) by 
the formula 

(70) logS{C) = - [ logmiga{p)-a{q))da{p), 

Jc 
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where q can be chosen to be any point in C. 

We will often need to consider expressions of the form ||t?||(5(a(p) — a{q)) 
for p, q E C, so to make things more readable we will from now on drop the 
a from the notation. 

Theorem 24 ([dJOU Theorem 2.1.2]). Letp, q ^ C withp not a Weierstrass 
point. Then 

(71) Gip.,y = sicy"' "!'''f 

where G = exp((7) is the exponential Green's function with respect to the 
canonical (1,1) form, and the product in the denominator is over the Weier- 
strass points of C taken with multiplicites, which are all g{g — l)/2 in the 
hyperelliptic case. 

The formula still holds for p a Weierstrass point; in this situation both 
the numerator and denominator vanish to order g[g — l)/2, and we must 
take the ratio of the leading coefficients of power series expansions about p 
in the numerator and denominator. 

10.3. Proof of Lemma 1201 To prove the Lemma, it suffices to show the 
following: 

Proposition 25. Given p > as in Definition \18l there exists a constant 
c{p) such that for all pairs of pairs {p,Vp), {q,Vq) with p G Vp, q G Vq, Vp, 
Vq elements of the cover V, we have 

(72) log|G(p,g)/G(po,go)| <c(p), 
where po, qo are the centres ofVp and Vq respectively. 

The proof of this proposition falls into two cases, depending on whether 
or not both pQ and qo are Weierstrass points. 

10.4. Case 1: not both Weierstrass points. Suppose at least one of 
Po and go is not a Weierstrass point. Since the Green's function is symmet- 
ric, we may assume without loss of generality that po is not a Weierstrass 
point, so both the numerator and denominator of the formula in [M] are non- 
zero. By continuity of theta functions, we may shrink p so that neither the 
numerator or denominator vanishes on the closed sets Vp and Vq. 

Now it is clear that the expression 



(73) 



G{p,q) 



G{po,qo) 



is bounded away from zero, and so it remains to give an algorithm to find 
such a bound. We do this by bounding the dervatives of the various quan- 
tities involved. 

10.4.1. Bounds on hyperelliptic integrals. In order to apply bounds on 
the derivatives of theta functions to bounding the above expression, it is 
necessary to understand the local behaviour of the map 

(74) a:C^ C/A 

given by hyperelliptic integration. For convenience, we give 
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(1) a choice of path from po to p, 
and 

(2) an exphcit bound on the resulting value |a(po) — ct(p)|. 
We briefly recall the setup of hyperelliptic integration: let ujq 



d X 

y ' 



9-1 



be a basis of differential 1-forms on C, let {Ai,Bi 



a symplectic homology basis, and let {wj 
normalised basis of differential forms such that 



1, 



l,...g} be 
. . ,5} be a 



(75) 



A, 



Let A C C^^^ be the resulting period matrix. Let S C be a fundamental 
domain, and let a : C(C) — )• S be the map obtained by integrating, aj will 
denote its ith component, obtained by integrating Wj, and similarly we set 



(76) 



ai{z) 



UJi 



It is possible to compute Qj, A and D and to evaluate a at given points 
to high precision due to work of Paul Van Wamelen implemented in MAGMA 
|B(]P97j . 

Proposition 26. For each disk V £ V, we can find an explicit compact 
box a[V] in D which contains the image a{V) of V under a, such that the 
diameters of such boxes tend to zero as p — ?• uniformly in V . 
Write p for such a uniform bound on the diameters. 

Proof. Since the are given by known linear combinations of the a-j, it 
suffices to show this result for the a.i. Let pq be the centre of V , and p 
the radius of its projection to P^(C). Without loss of generality, we assume 
I^Cpol — kpol- Throughout this proof, given X G C, Y{X) will denote a 
square-root of f{X), chosen to be continuous along radial paths if po is a 
Weierstrass point, and otherwise chosen to have no branch cuts in V (the 
cover V was carefully constructed so that this is possible). 

Case 1: po not a Weierstrass point 
Fix p £ V. We parametrise the path 7 = 7po,p by X{'y(t)) = + (Xp — 
Xp,))t. Thus for alH E {0, . . . , g — 1} we have 



\ai{po) - ai{p)\ < 



X' 



(77) 



< 



< 



Y{X) 
X^ 



dX 



Y{X) 



dX 



Y{x{jm 

Irl 



W{t)\dt 



< p sup 



Case 2: pQ is a Weierstrass point 
Fix p £ V. We parametrise the path 7 = 7po,p by X{'y{t)) = + (Xp 
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XpJt3/2_ rj.^^^ alH G {0, . . . , g - 1} we have 



\a{po)i - a{p)i\ < 



Y{X) 



dX 



(78) 



^ 2 I^Po ^pl 



< -/? sup 

2 reBp(Xj,„) 



YiXm)f 



Y^/(r)/(r -po) 



It is easy to check that the given bounds tend to zero uniformly with p, 
so we are done. □ 

10.4.2. Bounds on the derivatives of theta functions. To complete the proof 
of Case 1 of Proposition \25\ we need the following lemma on the derivative 
of theta functions: 

Lemma 27. Fix e > 0. Let z, w & such that z lies within distance e of 
the fundamental domain T), and w lies within distance e of z (both distances 
in the metric). Let c\ and C2 be positive constants such that 

ci 



(79) 



max y < 



27r' 



and 9(A) > C2. 



Set t{n) = IT Clin — ^^^), and write 

^ r 27rc2 1 



A = exp( 



47rc2 ' 



E 

n=0 



2^ ' 



and 



B = 27re'*''<=2 



V 



/vrc2 



n=0 



n)e 



1 ^ 



2JC2'K 27rC2 



r_£i_-| 

27rc2 I 

E 

n=0 



+ 



Cl 



47rc 



3/2 



two constants independent of e. Then we have the following hounds: 



1 



dzi 



< 23A3-'B 



2) \'&{z) -^{w)\<e23A3-^B. 
Suppose also that \'&{z)\ >c> l^A^-^B. Then 



3) 
4) 



d{l/d) 



dzi 
1 



1 



< (c - e23A9-^By^29A3-^B 
< e(c - e2<^A3-^B)~^29A9-^B. 



Proof. Detailed background material for this proof may be found in [Mum83[ 
II §1]. From the power series expansion 

t?(z) = exp {mn^Kn + 2'Kin ■ z) , 



(80) 
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(here i denotes a square-root of —1, to distinguish it from the index i) we 
see 



9ir 



(exp (TTin^An + 2inn ■ z) {2mnj)' 

neZs \ V i i / J 

< 2^ X I cxp ( -7rc2 X + ci J (27rn,)'- j 

— ( X] ^■'^P (— 7rc2n^ + cin) J I y~^(27rn)'' exp (— 7rc2n^ + cin) 

VneN / VneN 



Now 
(81) 



x^e ^ dx = -r 

2 V 2 



so recalhng t{n) = ■^jTTCQ,{n — 2^)) we obtain 
exp(— 7rc2n^ + cin) 

nGN 



<exp V exp(-t(n)2)+ / exp(-i(n)2)d 



n 



' I 27rc2 



r_EL-| 

I 27rc2 ' 



exp 



, 2 \ / ^'"^2 ' 



n=0 



2^^ 



Now wc must do the same for ^^gpj(27rn) exp(— 7rc2n^ + c\n) (exactly the 
same argument would work for the rth derivative, but we need only the first 
derivative) : 

(27rn) exp(— 7rc2ra^ -\- c\n) = 27rexp( ^ ^ ) n exp(— t(n)^) 



neN 



(ii) (vfe ^ ^^"^ ^"^^-^^"^'^ + ^ ^ exp(-.(n)^)] 
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Now 
(82) 



^t(n)exp(-t(n)2) < ^ t{n) exp{-t{nf) + \t{n)\ exp{-t{nf) d 

< V t(n) eyip(-t(nf) + ^= texp(-t'^)dt 

±-1 xAc2 7t=0 



n 



n=0 



t(n)exp(-t(n)2) + 



n=0 

SO combining the above results we find 
^^(27rn) exp(— 7rc2n^ + cin) < 

neN 

2vrexp(-^ 

V47rC2/ y/TTC2 



2-v/vrci' 



1 



I 27rco I 



-i(n)2 



V 



n=0 



2v/c27r 27rc2 



+ 



Cl 



n=0 



47rc 



3/2 



This concludes the computation of the bounds for (1). The bound for (2) 
is simple; it is simply e times the bound for (1). For (3), we use 



(83) 



a(iA?) 

dzi 



(^) 



-1 dd 



(^) 



i9(z)2 dZi 

to conclude that the bound for (3) is (c — e2^A^^^i3)~2 times the bound for 
(1), and similarly that the bound for (4) is e times the bound for (3). 

□ 

10.5. Case 2: po and go both Weierstrass points. In this case, both 
the numerator and denominator of (j7ip vanish to order g{g — l)/2 at pq. 
To approximate the value of the expression at (poiQo) is not hard using a 
limiting approach, but to control the derivative near to that point would be 
more difficult. Instead, we apply the maximim modulus principle and take 
advantage of the symmetry of the situation to reduce to Case 1. 

As before, write Vp^ and Vq^ for the non-overlapping disks centered at 
Po and go in V. Using the maps induced in Section 19.2.21 we shall abuse 
notation by viewing Vp^ and VJ^q as open discs in C. Recall that a : C — )• 2) 
is the map to a fundamental domain inside obtained by integrating. Then 
write 

C 



• ^0 ^ ^qo 



(84) 

and 
(85) 

K{p, 



{p, q) ^ 



i){ga{p) - a{q)) 



5(X)det(9(T 



exp(-7r9(ga(p) - a(g))*(S'(r))~^Im(ga(p) - a{q))) 
Uwew exp(-7r3=(c/a(p) - a(u'))*(9(r))-i/m(ga(p) - a{w))) 
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where the product is over Weierstrass points w counted with their weights 
(which we recall are all g{g — l)/2 in the hyperelliptic case). 
Then (p is related to the gth power of the Green's function by 

(86) G{p,qy = \4>{p,q)\K{p,q), 
and so to bound 

(87) Gip,qy, 

it suffices to bound \<j){p,q)\ and K{p,q) separately. 

Bounding K{p, q) is straightforward; Since ^(t) is positive definite, 9(r)~^ 
is also, and so it is sufficient to ensure that the values of ga{p) — a{q) and 
ga{p) — a{w) do not get too small. This is achieved by making a suitabel 
choice of fundamental domain D. 

To bound 0(p, q) on V^^ x Vgg is a little harder, since it is written as a 
quotient of two holomorphic functions which both vanish to the same degree; 
bounding the numerator and denominator separately will not be sufficient. 
Instead, we will apply the maximum modulus principle, using the following 
lemma: 

Lemma 28. Retaining the setup above, there exists a computable union of 
eight subsets Vj C C x C with the following properties: 

1) sup(p_^)gy^^xy^^ |(/>(p, q)\ < maxig|i^,„^8| sup(p^g)gv, <P{P - Q)- 

2) For each i, at least one of p and q is never a Weierstrass point for 

(p,g) e V,. 

Proof. The construction of the Vj is simple; they are the eight 3-real-dimensional 
faces of a 4-real-dimensional cube contained in C x C, centered at {po, go) and 
chosen just large enough to contain V^^ x Vg^^. If is easy but un- illuminating 
to write down parametrisations of the faces. 

Assertion (1) now follows immediately from the maximum modulus prin- 
ciple applied to this cube. 

It is easy to check that this closed cube does not contain any images of 
Weierstrass points apart from po and qo, and on each face at least one of the 
real or imaginary parts of either p or g is constrained to be unequal to the 
corresponding coordinate of po or qo respectively, so (2) follows also. 

□ 

Using this lemma, we see that we may apply the derivative bounds from 
Lemma [23 to enable us to numerically bound the supremum of \4>{p,q)\ for 
{p,q) G Vj; on a fixed Vj, we know from (2) that at least one of j? or g 
is never a Weierstrass point, so (posibly after interchanging p and q), the 
expression for cp is given as a ratio of non-vanishing holomorphic functions, 
whose derivatives are bounded in Lemma [271 and which we may therefore 
rigorously numerically bound. 

This gives us an upper bound on \(j){p,q)\ for (p, g) G Vpf^ x Vqg. We also 
need a positive lower bound, but this is easily obtained by applying the same 
algorithm to the function l/0(p, g). 

11. The first naive height 

Recall that K = Q. 
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Definition 29. We begin by defining a naive height H : A[K) — >■ K>i. 
Note that it does not arise directly from a projective embedding, but we still 
call it 'naive' since it is relatively simple to define and compute, and is an 
approximation of the Neron-Tate height (as we will shortly prove). 

Let iJ,> be less than half of the shortest distance between any two Weier- 
strass points under the Archimedean absolute value of K . 

Given p G A{K), write p = [D — deg(D)oo] where D is a reduced divisor 
on C . If the support of D contains any Weierstrass points, replace D by the 
divisor obtained by subtracting them off; this equates to translating p by a 
2-torsion point, and so will not affect the Neron-Tate height. Let d denote 
the degree of the resulting divisor D . 

Choose once and for all a pair of degree-d effective divisors oOp and oOp 
with disjoint support, supported on Weierstrass points away from oo, such 
that no point in the support of D is within Archimedean distance ^ of any 
point in the support of oOp or oo^. The existence of such divisors is clear 
since there are 2g-\-l Weierstrass points away from oo and reduced divisors 
have degree g. 

Recalling that D~ denotes the image of D under the hyperelliptic involu- 
tion, define 



where L/K is the minimal finite extension over which D, oo^ and oOp are 
pointwise rational, recalling that if D = di, oOp = Y^- qj and oOp = qf 
then 

/ooA ^ /n In- 2^ Yj (UjH^Pj)du{qhq^ 

(89) d„(D-oo„Z) -°°.)-n ,^(,,,,.),^(,-,,;) . 

We define a logarithmic naive height by J^{p) = \og{H(j))). 

Remark 30. Note that in the above definition, we made use of the assump- 
tion that K = Q - indeed, this is the only point at which we use this fact. 
The condition would not be hard to remove by allowing different choises 
of divisors oOp and oOp at different Archimedean absolute values, and then 
showing that the change to the resulting height was not too large. 

Proposition 31. The products in the definition above are finite; in partic- 
ular, the heights are well defined. 

Proof. Prom the definitions of the metrics over non- Archimedean absolute 
values, it is clear that di,{D — oOp,D~ — oOp) = 1 for all but finitely many 
such absolute values. □ 

Combining previous results, we obtain the following theorem, which is the 
main result of this paper. 

Theorem 32. For all p G A{K) we have 



(90) h{p)-J^fip) 



< ^1 + ^2 + . 
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where is from Definition ^2 is from Corollary [3 and =^3 is from 
Definition [21[ 

Write c for the constant SS\ + SB2 + <^3 • 

12. Refined naive heights 

We define two new naive heiglits whicli are eacli in turn simpler to com- 
pute, and we bound their difference from the Neron-Tate height. The last of 
these heights will be sufficiently simple as to allow us to give an algorithm 
to solve the saturation problem. 

Definition 33. Given p G A{K), we take the divisor D = Yli=iPi over 
some finite L/K as in Definition{2^ and write pi = (Xj, Yi). Then set 

d 

(91) \P{p) = Y,h{X,^), 

i=l 

and set 

(92) ht(p)=h(^n(X-X,J^, 

where the right hand side is the height of a polynomial, which equals the 
height of the point in projective space whose coordinates are given by its 
coefficients. Given B > 0, define 

(93) M'^{B) = {peA{K):h'^{p)<B} 
and 

(94) (B) = {p e A{K) : h\p) < B} 

Our main aim in this section is to give computable bounds on the differ- 
ences 1=^ — h'^l and |h^— h^|. 

Lemma 34. There exist computable constants < ci < C2 with the following 
property: 

for all non- Weistrass points p = {x : s : y) £ C{K), and for all Archimedean 
norms \ — on K, we have 

(95) ci < d,(p,p-)/(2min(|y|^ , \Y'1)) < C2, 

where as usual we write Y = y/s^^^ and Y' = y/x^^^. 

Proof. Fix an Archimedean absolute value v. Recall that d = d;^ is the 
weak-pseudo-metric given in Definition [T2] using a Merkl atlas. There are 
two cases to consider. 

Case 1: there is no chart U in the Merkl atlas such that p G [/''^ and 
p~ G C/^i. Then the value of d{p,p~) must be one of a finite set of positive 
values (as defined in [TS|) . We also note that min(|y|j^ , |y |^) is bounded 
above (for obvious reasons) and is also bounded below since p cannot be too 
close to a Weierstrass point (otherwise p and p^ would both be contained in 
the same Merkl chart), and so min(|y|^ , |i^'|j,) is bounded away from zero, 
to the result follows. 

Case 2: there is a Merkl chart U with p G [/''^ and p^ G C/'"^. Clearly 
such a chart must contain (and so be centred at) a Weierstrass point, call it 
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w = (xu) '■ Su, ■ Uw)- Without loss of generality say \xyj\ /le so there is a 
t>l such that \xp/sp\ < t. Then |F| < so it suffices to show that 

|y| is close to d{p,p~). Now d{p,p~) is given by the following procedure: 

Translate w to the point X = using a Mobius transformation by an 
orthogonal matrix 



(96) 

Then d{p,p-) = 2 
(97) 

whereas 
(98) 



M 



— Xy^ 
-1 



. Thus 



dip,p-) = 2 



Xp Xyj 



2|Fp| = 2 



where fw{t) = f{t)/{t — Xy^). Thus by bounding the rational function 

(99) {X - Xy,)/U 

on the domain (recalling that R is the constant defined in 4.2.1) 

(100) {X eC:\M -XlKR} 

(a region on which, by construction, it has neither zeros nor poles), we are 
done. □ 

Definition 35. Let L/K he a finite extension, and let p ^ q E C{L) be 
distinct points. Set 



(101) 



(P, 9)l = 7^7^ log Yl Mp^Q)- 



Lemma 36. There exists a computable constant c with the following prop- 
erty: 

let L/K be a finite extension, and let p = {x : s : y) E C{L) be a non- 
Weierstrass point. Then 

(102) \{p,p-)^-ig + l)hix/s)\<c 

Proof. For | — | non-Archimedean, wc have that if |a;| < \s\ then d{p,p~) = 



\2y/s^^^\, and if |s| < |x| then d{p,p ) 
Archimedean we obtain 



|2y/x^"'"^|. Hence for non- 



|2t/|^min(l/ \x\l'^^ , 1/ \s\ 



(103) d^{p,p , 

We have shown above that for Archimedean v we have computable < ci < 
02 such that 

(104) ci < d,{p,p-)/mm{\2y/x3+^ , \2y/sa+^) < C2. 

Hence 
(105) 



n vc2< 



< 



n vci. 
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Now rii^eMf" ^1 '^^^^'^^ is bounded uniformly in L, and similarly for C2. Fi- 
nally, note 



(106) 



□ 



Definition 37. Assume that the hyperelliptic polynomial f is manic. Given 
a Weierstrass point d with ^ 0, set fd to be the univariate polynomial 
such that for all p ^ d E C{K), we have 

(107) fd{Xp){Xp-Xd) = f{Xp). 

It is clear that fd will have integral coefficients, since f does and Xd is a 
root of f. 

Lemma 38. Let L/K be a finite extension, and let p, d E C{L) such that 
Sp ^ and d is a Weierstrass point with Sd ^ 0- Assume further that 
the hyperelliptic polynomial f is monic, so that Xd is integral. Let v be a 

non- Archimedean absolute value of L, and suppose \Xp — Xd\^ < 

Then 

|2 



fd{Xd] 



(108) 



\Xp — Xd\ 



fd{Xd 



Proof. By definition, we have 



(109) 



= - ^d\t, fd{Xp) 



SO it suffices to show that 



fd{Xp) — fd{Xd 



. Writing 



(110) fd{Xp) = {Xp - XdT + ^{Xp - XdT-' + • • • + *(Xp - Xd) + fd{Xd) 

where the coefficients ★ are integral, we sec that the greatest norm of any 
term on the right hand side is achieved by fd{Xd) and no other term, so the 
result follows. □ 

Lemma 39. Let L/K be a finite extension, and let p ^ d e C{L) be such 

that Sp ^ Q and d is a Weierstrass point with 7^ 0. Assume further that 
the hyperelliptic equation f is monic, so that Xd is integral. Then 

(111) - \ogd,(j>,d)<[L:K](^]i{Xp-Xd)+W{Xd))\. 



Note that the sum is over the non- Archimedean absolute values. 

Proof. The right hand side naturally decomposes as 
(112) 

[L:K]C-\i{Xp-Xd)+hCfd{Xd)))= \log'\Xp-Xd\-^+\og+\fd{Xd) 
Now it is clear that 



(113) 



Y \\og+\Xp-Xd\-;^+\og+\fd{Xd) 



>o. 
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SO it suffices to prove that for each non- Archimedean v we have 



(114) 



log(d,(p, d)) < i log+ |Xp - Xa\-^ + log+ |/rf(Xrf) 



or equivalently that (at this point we drop the subscript v from the norm) 



(115) 



-1 



d,(p, d)-' < max{\Xp - X^r' , 1) max( UXa) , 1)' 



Recalhng that 



fd{Xd] 



for simphcity, we see 



< 1 and writing F 

this is equivalent to showing 

(116) d,(p,d)2 >F2min(|Xp-Xd|,l). 

We divide proving this in to two cases. The first is when \Xp — Xd\ > F. 
Then 

(117) 

so Equation (jll6|) fohows. 

The harder case is when \Xp — < F. We apply Lemma [3S] to see that 
\Yp\' = \Xp - Xd\ F, and so 



(118) 



> 



m8.x{\Xp - Xd\\\Y/) if \Xp\ < 1 
1 if \Xp\ > 1 

ma.x{\Xp-Xd\\\Xp-Xd\F) if < 1 
1 if \Xp\ > 1 

Fmax(|Xp - Xd\^ , \Xp - Xd\) if \Xp\ < 1 



1 



if \Xp\ > 1 



F \Xp - Xd\ if \Xp\ < 1 
1 if \Xp\ > 1 

>Fmm{\Xp-Xd\,l) 

>F^mmi\Xp-Xd\,l) 

Lemma 40. Let Xi, X2 G L. Then 

(119) H(Xi + X2) < 2*^^lI{L:K] Yi^Xi) H(X2). 

Proof. 
(120) 

H(xi + X2)[^^^i = n + 

< I n max(l,|Xi|,,|X2|J I I n max(l,|Xi|, + |X2|J 



n 



< H max(l,|Xi|Jmax(l,|X2|J J] 2max(l, max(l, IX2IJ 



. ueM? 
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as required. □ 

Lemma 41. Fix fi > and n S N. There exists a computable constant 
4>^{n) with the following property: 

Let L/K he a finite extension of degree at most n, and let p, d £ C{L) 
such that Sp ^ and d is a Weierstrass point with Sd ^ 0. Assume further 
that the hyperelliptic equation f is monic, and also that for all Archimedean 
absolute values v E M'^ , we have du{p,d) > /i. Then 

(121) {p,d)L<^h{Xp) + 4>^{n). 

Proof. Prom Lemma [39] and the fact that di, {p,d) > fi for all Archimedean 
v, we see that 

(122) {p, d)^ < \ h(Xp - Xd) + h(/,(Xd))) - log(/.)#Mr /[L : K\. 
Now by Lemma HOt we have 

(123) h(Xp - Xd) < HXp) + h{Xd) + log(2) , 
so for fixed L and d we may take 

(124) 4>t^''{n) = - HfdiXd))) - log(/i)^^ + h(X,) + ^^-^ log(2). 



Thus the existence of a bound uniform in L and d is clear (since there are 
only finitely many Weierstrass points, and the degree of L is bounded). □ 

Lemma 42. There exists a computable constant c such that the following 
holds: 

given p G A{K), let D, oOp and oOp denote the divisors given in Definition 
Let L/K be the minimal finite extension such that D, oOp and oOp are 
all pointwise rational over L. Then we may write 



d 

D = Y.P. 



i=l 
d 

(125) ^l = Y.1^ 



i=l 
d 

2 



OOp 



i=l 



Then 



(126) ^(p) > ^ {p.,pr)^ - ^ {p,,qj)^ - ^ {p„q'^)^ + c. 

1=1 \ j=i j=i J 

Proof. Recall that 
(127) 

d d d d 

i,j=l i,j=l i,j=l i,j=l 
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Since the qi and q'^ are distinct Weierstrass points we easily bound J2ij=i (^Qi^ ^ 
For i ^ we have by Theorem [T6l that 

(128) (P^'PJ)^ ^ -^(M) • #Afr/[L : K] - ^i, 

where M is the Merkl atlas constructed in Section 19.21 (since all other con- 
tributions to the pairing are non- negative), so the result follows. □ 

Lemma 43. There exists a computable constant d such that in the setup of 
Lemma\4^we have 

d 

(129) J^{p)>^h{Xp^) + c' 

i=l 

Proof. In Lemma H2] we showed 

d / d d \ 

(130) jr(p) > ^ {m,P7)L - E (P^^ ii)L - E {p^^i'^)l + ^• 

i=i \ j=i j=i j 

In Lemma [36] we showed (using that the pi are never Weierstrass points) 
that for some computable c\ we have 

(131) |(K,pr)^-(5 + i)h(Xpj| <ci. 

Write n for a positive integer such that the minimal field L appearing in 
Definition [29] can always be taken to have degree at most n (this n will 
depend only on the genus of the curve C). In Lemma [H] we showed (using 
that d;y(pj, qj) > n where /i is as in Definition [29]) that 

(132) {pi,q,)^<^h{XpJ + <p^{n), 

and similarly for q'j. 

Combining these, we see using d < g that for each i 

(133) 



da 1 

E (P^^ - E iP^^ '^ia + l) HXn) - 2 E 2 ^^^^'^^ " + ^'^'^'^ 
j=i j=i j=i 

= (((7 + 1) - 2di) h(XpJ - ci + 2#^(n) 

> h(XpJ - ci + 2d,^^(n). 

from which the result follows. □ 

Theorem 44. There exists a computable constant c such that for all p G 
A(K) we have 

(134) h{p)+c>li^{p). 

Proof. From Theorem [H2] we know that there exists a computable constant 
c' such that 

(135) h{p) + c' >J^{p). 

The result follows from combining this with Lemma [l3l □ 
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Corollary 45. For any constant B: 

(136) M{B) C M'^iB + c) 

where c is the computable constant from Theorem 

Lemma 46. Fix a finite extension L/K. Given ai, . . . , o„ G L, set ip„ 
n"=i(*-«t)- 'Then 



(137) 



i=l 



< #M^log(4)(n2+n-2)/2. 



Proof. From |Lan831 Chapter 3, Proposition 2.4] we have for ah m > 2 that 

(138) |h(i - am) + Hi^m-i) - h(V'm)| < m#M^ log(4) 

(note the difference in normahsations between our heights and Lang's). The 
formula follows by induction and using that h(i — Oj) = h(aj). □ 

Corollary 47. For all p G A{K) we have 

(139) |h^(p)-ht(p)| <#M^log(4)(5'+g-2)/2. 
The main result of this chapter is then 



Theorem 48. Let c be the computable constant from Theorem 44- Then 
for all constants B we have 

(140) M{B) C Aft (5 + c + log(4)((72 +g- 2)/2) . 

The point is that these finite sets M^{B) are effectively computable, so we 
can in turn use the results from |Holl2aj to compute the finite sets M{B). 
We describe one algorithm to compute M\B), setting K = for simplicity: 

1) Let S be the finite set of all polynomials nf=i(^ ~ d < g, of 
height up to S. 

2) for each polynomial a G 5, if a is not irreducible remove it from S and 
insert into S each of the irreducible factors of a. 

3) it suffices to determine for each a £ S whether a is the 'x-coordinate 
polynomial' of a divisor in Mumford representation; in other words, whether 
there exists another univariate polynomial b such that (a, b) satisfy the prop- 
erties of a Mumford representation. Now the polynomial a also determines 
a set of 2deg(a) distinct complex points on C - the preimages of zeros of 
a under the hyperelliptic projection. These can be computed to any finite 
precision. Such points will satisfy y = b{x), thus if we can bound the de- 
nominators of the coefficients of b then we can find a finite precision to 
which we need to compute the complex points to see if they correspond to a 
polynomial b with rational coefficients. Such a bound on the denominators 
is supplied by the following proposition. 

Proposition 49. Let K/Q be a finite extension, with integers Ok and p a 
prime ideal in Ok- Let g > be an integer. Let f,h£ K[x] be polyno- 
mials which are integral with respect to p (ie the p-adic valuation of their 
coefficients are positive) and such that 

- Af + h'^ is separable; 

- f has degree 2g + 1; 
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- h has degree at most g + 1. 

Fix an integer 1 < d < g. Suppose we are given a pair of polynomials 
a = Yli=o ^i^* of degree d and b = YliZo ^i^* degree at most d — 1 in K[x] 
and a constant c € K with the following properties: 

- a, b and c are integral at p (ie their p-adic valuations are positive); 

- a is primitive; 

- ordp c > mini(ordp 

- A = disc(a) is non-zero. 
Suppose also that 



(141) 
Then 



(142) ordp c < - ordp A + 



b 

c I \ c 



h-f. 



d — d + max 



2g+l 



, deg{h) ordp ad- 



Proof. Let L be a 'sufficiently large' extension of the completion Kp] by this 
we mean that L is a finite extension of the completion which we will require 
to be closed under taking roots of a certain finite collection of polynomials, 
which will be described as we go along. Let vr denote a uniformiser of L, 
and Ol its integers. 

We assume a splits in L; write xi, . . .Xd for the roots. Each Xi may be 
uniquely written j/vr^' where rj > and Xj E Ol has minimal 

valuation. Let 



(143) 



a = JJ(7r'''x - Xi) eOiix]. 



i=l 



Now a and d have the same degree and (distinct) roots, are integral, and 
both have at least one coefficient which is a unit in Ol, hence a and d differ 
by a unit in Ol. 
Let 









l}ri 










1)t2 
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M = 














l)r-d 




so det(M)2 


= disc(a) 


= unit X A. 


Let 








/ 1 


Xi 


(145) 






1 


X2 




M = 










I 1 





'X2 



'Xd 



Xi 
X2 



By ()14ip we know that for all 1 < i < d we have 

d-i 

(146) 

j=0 



E (7) = 



xr 



d-1 



d-1 
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for some yi in L (obtained by assuming L sufficiently large) satisfying yf + 
h{xi)yi = f{xi), and hence that 



(147) 



/ ^0 \ 



y2 



\yd J 



In order to bound above the order of c at vr, it therefore suffices to bound 
below the order of the right hand side of (jl47p at vr. We do this in two steps. 
Firstly, we easily obtain from properties of valuations that 



ord,ryj > -r^max ( deg(/i) ) . 



(148) 



Secondly, we must do the same for . Now M ■ Ti-Ccf-i) max; ^ matrix 
over Ol, and hence so is its transposed matrix of cofactors, which we shall 
denote Mc- We then find that 

(149) "^-1 ^ ^ 



M 



TT' 



d{d—l) maxi 



det(M) 



X M, . 



We also compute that 



(150) 

and hence 
(151) 



det(M) = ^-(''-1) S'^'' det(M) 

= 7r-('^-i)S'^VA-unit, 



M- 



a/A • unitvr'^^'^"^) max,(ri)-(d-i) '£^ n ' 



We also note that = ord,r cLd, and maxj(rj) < J2i so combining the 

above results we see that 

(152) ord^ ~ 2 ^ ~'~ i^'^ ~ ^^"^ ^ ^scx. ^ ^ ~^ — , deg(/i)^ ^ ordjr a^, 
from which the result immediately follows. □ 
We note that in the elliptic case, we recover the classical result that \ a\. 



13. Appendix: an algorithm to gompute a gover of the Riemann 

SPHERE 

In this appendix, we will describe an algorithm to determine a cover of 
the Riemann sphere by open discs in the chordal metric, satisfying certain 
technical hypotheses coming from Merkl's theorem. 

Let S"^ denote the Riemann sphere, and fix r < 1. By an 'algebraic point 
of 5"^' we mean a point which maps to an element of Q'^^s C. Note that 
such a point can be described by a finite amount of data. 

Definition 50. We say a closed subset Z C S'^ is 'starry' if the complement 
of Z is a finite union of open discs (in the chordal metric) of radius r < 1 
centred at algebraic points. 



36 



DAVID HOLMES 



Proposition 51. Suppose Z is starry, has no isolated points and the interior 
Z° is non-empty. Then there exists an algorithm to find an algebraic point 
p E Z° such that Z\Br (p) has no isolated points. 

Proof. The set of points p such that Z\Br{p) has isolated points is contained 
in the set of points p such that there exist qi, q2 centres of discs in the 
complement Z'^ and t e Z with 

(153) d{p,t)=d{t,qi)=dit,q2)=r. 

The set of such points p is contained in a finite union of computable circles 
(boundaries of discs) of radius r and centred at algebraic points, and so the 
complement of that set in Z° is non-empty open and described by a finite 
collection of computable conditions, and we can find an algebraic point inside 
it. □ 

Proposition 52. Let Z be starry. Then every connected component of Z 
has either non-empty interior, or is an isolated point. 

Proof. Fix a connected component Y. If Y° = 0, then Y = dY, so y is a 
finite union of circles of radius r. Since the complement of Z can be covered 
by a finite collection of open discs Di of radius r, we see that there must be 
a segment of Y of positive length which is touched by exactly two of the D^. 
These discs therefore have the intersection of their boundaries containing 
a segment of Y of positive length, which is impossible since the discs are 
distinct and have radius r < 1. Compactness of Z then shows that Y is 
isolated. □ 

Definition 53. Given a disc D and an integer n > 3, let denote a 

smallest (unique up to rotation) regular n-sided polygon (with sides consist- 
ing of arcs of great circles) containing D, and similarly the largest such 
polygon contained in D. Note that if D has algebraic centre, then and 
-D"" have algebraic vertices. 

Given a finite set of algebraic points S, set 

(154) U{S) = {Br(jp):peB}, 

(155) Un{S) = {{Br{p))n : P & B}, 

(156) U^iS) = {{Brip)r:peB}. 

Lemma 54. Given a finite set S of algebraic points, the following hold: 

1) ifU{S) is not a cover, and the complement ofU{S) has no isolated 
points, then there exists n > such that U"'{S) is not a cover. 

2) if U{S) is a cover, then there exists n > such that Un{S) is a cover. 

Proof. 1) If U{S) misses a point, then it misses an open set, in particular a 
disc of some radius e > 0. Eventually, the sets for D G U{S) will have 
diameter less that 2r + e/4, and so {/""{S) will miss an open set of points. 

2) this is clear by a similar argument, once we observe that given any 
cover of a compact set X by open sets {Ui}i, there exists some e > such 
that for every point p e X there exists an index i such that 

(157) B,{p) c Ui. 

□ 
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Proposition 55. Given a finite set S of algebraic points and an integer 
n > 1, there exists an algorithm to determine whether U"'{S) is a cover, and 
to determine whether Un{S) is a cover. 

Proof. It suffices to give an algorithm which, when given a finite set of 
polygons (with arcs of great circles as their sides), will determine whether 
these cover 5*^. This is easy; one approach is to subdivide the union of the 
polygons into a finite set of triangles which meet properly, then to compute 
the homology of the resulting complex; we have a cover if and only if the 
second homology group is non-zero. □ 

Proposition 56. Fix B a finite set of points in 5^, and choose a rational 
number < r < 1 less that half of the shortest distance between two distinct 
points in B. Then Algorithm\57\ will yield a finite cover C of S'^ by discs of 
radius r satisfying the following properties: 

a) B is contained in the set of centres of discs in C; 

b ) Given a disc D in C and a point p £ D which is the centre of a disc 
in C, then p must be the centre of D; 

c) there exists a computable e > such that for all p, q centres of discs in 
C, d{p,q) > r + e (equivalently, property (1) still holds when the discs are 
replaced by their closures). 

Algorithm 57. 1) Set S = B; 

2) Set N = 1, and mi = 1; 

3) forl<n< N: 

if Um„ (S) is a cover, output S, finish; 

else if [7™" (5) is not a cover, use Proposition [57] to find an algebraic point 
P{N,n) outside U^"{S) such that the complement of U{S L) {p(^]\f n)}) ^'^^ 
isolated points; 
else mn+ := 1; 

4) N+:= 1, ruN := 1. Go to (2); 

Proof. It is not hard to see that if the procedure terminates then the resulting 
cover has the desired properties, but it remains to prove termination. We 
argue this by contradiction; suppose that the procedure does not terminate. 
This yields a cover of S"^ by U{S) for some infinite set of algebraic points 
5, so by compactness there exists a finite subset T C S such that U (T) is a 
cover. Let Nq > be such that 

(158) TCBU {p(,,„) : / < n < iVo} 

By Lemma [SH there exists no > 1 such that C/„g(T) is a cover, and so the 
algorithm will terminate when N = Nq + uq. □ 
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